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1.  Introduction 

The  dynamic  effects  of  combustion  arise  from  the  thermomechanical  response  of  a  com¬ 
pressible  medium  to  localized  heat  addition.  Thermal  expansion  of  a  fluid  element,  caused 
by  a  rapidly  rising  temperature,  induces  wave-phenomena  into  the  surrounding  material.  One 
may  observe  acoustic  fields,  weak  or  strong  shock  waves,  depending  upon  the  time  scale  of 
heat  addition.  The  interaction  between  the  chemical  processes  and  the  induced  gasdynamic 
phenomena  are  essential  to  an  understanding  of  ignition,  explosion,  high  speed  deflagrations, 
detonations,  and  transition  of  flames  to  detonations.  These  topics  are  integral  to  the  study  of 
accidental  explosions,  premature  ignition  of  explosives,  blast  damage  caused  by  large  fires  and 
internal  combustion  engine  knock,  to  name  a  few  examples. 

During  the  3-year  period  of  our  research  program,  dealing  with  the  dynamics  of  combus¬ 
tion,  we  have  focused  on  the  following  areas: 

(1)  Source  theory  for  shock  wave  generation  in  inert  and  reactive  compressible  gas  mix¬ 
tures  including  detonation  wave  initiation. 

(2)  Evolution  of  an  autoignited  thermal  explosion  in  a  confined  reactive,  compressible  gas 
beyond  the  induction  period. 

(3)  High  speed  steady  combustion  waves  with  complex  chemical  kinetics. 

The  studies  in  (1)  and  (2)  are  fundamentally  transient  in  character  and  closely  related 
in  terms  of  conceptual  development.  They  comprise  an  organized  effort  to  model  the  basic 
physical  processes  occurring  in  dynamical  combustion  events.  Initial-boundary  value  problems 
are  formulated  for  each  system.  Solutions  describe  the  transient  response  of  the  system,  and 
permit  one  to  predict  which  type  of  combustion  process  evolves  on  relatively  long  time  scales. 
For  example,  interest  is  focused  on  the  factors  that  lead  to  flame  initiation  accompanied  by 
acoustics  in  a  given  configuration,  in  contrast  to  those  associated  with  the  generation  of  a 
detonation  wave  in  the  same  system. 

The  study  of  high  speed  combustion  waves  in  (3)  emphasizes  the  interaction  of  complex 
chemical  kinetics  with  a  steady  compressible  flow.  Chain  branching  and  termination,  as  well 
as  high  temperature  dissociation  and  recombination  are  incorporated  in  the  modeling  in  order 
to  emulate  processes  in  //j-Oj  and  hydrocarbon  systems. 

Mathematical  modeling  is  based  on  the  general  describing  equations  for  a  reactive,  com- 


2 


pressible  mixture  of  perfect  gases  including  transport  effects.  Formal  methods  of  analysis  are 
used  to  study  solution  properties  including  uniqueness,  long-time  evolution,  sensitivity  to  es¬ 
sential  parameters  and  bounds.  Limit  processes,  derived  from  perturbation  analysis,  are  used 
to  construct  rational  approximations  to  the  complete  equations,  for  various  parameter  values. 
Solutions,  constructed  in  terms  of  asymptotic  expansions  appropriate  for  problems  exhibit¬ 
ing  disparate  time  and  length  scales,  describe  the  complete  time-history  of  spatially  distributed 
processes.  This  approach  to  dynamical  combustion  problems  complements  the  complete  compu¬ 
tational  treatment  of  the  full  equations.  It  provides  general  parameter  dependence  of  solutions 
and  “benchmarks”  to  be  used  for  testing  complex  code  reliability.  The  latter  issue  is  of  im¬ 
portance  because  the  discretized  form  of  the  mixed  parabolic-hyperbolic  systems,  describing 
dynamical  combustion,  can  generate  a  variety  of  extraneous  nonphysical  phenomena  which  are 
difficult  to  detect  solely  by  numerical  testing. 

Progress  in  each  of  the  three  specified  research  areas  is  described  below.  More  detailed 
information  is  available  in  the  summary  descriptions  given  in  Section  2  or  in  the  publications 
themselves  listed  in  Section  3. 

Major  accomplishments  of  our  detonation  wave  initiation  work  include: 

(a)  development  of  a  robust  implicit  Navier-Stokes  solver  capable  of  providing  reliable  solutions 
for  reactive,  compressible  flows  with  transport  effects. 

(b)  characterization  of  initiation  of  detonation  in  a  reactive  gas  by  boundary  power  deposition. 

(c)  development  of  a  robust  MacCormak  scheme  for  the  reactive  Euler  equations  and  solution 
generation  for  volumetric  power  deposition  initiation  of  reactive  gasdynamics. 

(d)  in  preparation  for  studying  reactive  wave  initiation  in  compressing  gases,  mathematical 
methods  have  been  developed  to  study  the  dynamic  compression  of  a  gas  in  a  cylinder. 
Major  accomplishments  in  our  work  on  thermal  explosions  include: 

(a)  a  unified  theoretical  formulation  for  diffusive  and  nondiffusive  thermal  explosions. 

(b)  an  exact  solution  for  a  class  of  nondiffusive  thermal  explosion  IBV  problems,  which  inc- 
cludes  examples  of  global  rather  than  pointwise  blow-up. 

(c)  a  mathematical  model  for  reaction  transients  in  a  gas  with  dissociation-recombination 
chemical  kinetics. 

(d)  a  precise  description  of  where  thermal  explosions  are  located  for  gaseous  reactive-diffusive 
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ignition  models  showing  that  global  blow-up  can  occur  for  certain  types  of  reactions. 

(e)  a  description  of  how  the  thermal  explosion  singularity  evolves  for  rigid  reactive-diffusive 
models. 

Major  accomplishments  in  the  study  of  high  speed  steady  combustion  waves  include: 

(a)  the  modeling  of  initiation  of  hig  speed  reaction  waves  in  a  compressible  reactive  gas  with 
transport  effects. 

(b)  the  discovery  that  non-Chapman-Jouguet  conditions  exist  in  high  speed  reaction  zones 
when  the  chemistry  leads  to  reductions  in  the  mixture  mole  number  as  the  reaction  evolves. 


2.  Research  Activities:  A  Summary 

During  the  34-month  lifetime  of  the  research  project  described  in  this  Final  Report  the 
participants  published  15  articles  in  peer-judged  journals  covering  the  fields  of  combustion 
science,  fluid  mechanics  and  applied  mathematics.  In  addition  two  additional  manuscripts  have 
been  submitted  for  publication.  The  principle  investigators  and  their  students  have  given  a 
total  of  nearly  50  technical  presentations  at  colloquia,  meetings,  seminars  and  short  courses. 
Both  Profs.  Bebernes  and  Kassoy  held  NATO  Cooperative  Research  Grants  during  the  grant 
period  that  fostered  international  cooperation  with  research  colleagues  abroad. 

Outreach  activites  of  Prof.  J.  Bebernes  include: 

1.  Invited  talk,  EQUADIFF  6,  Czechoslovak  Academy  of  Sciences,  Brno,  August  2985. 

2.  Invited  talk,  Instituto  di  Matematica  Applicata,  Universita  di  Padova,  Padova,  January 
1986. 

3.  Visiting  Professor,  University  of  Wyoming,  January-June  1986. 

4.  Invited  Lecturer,  3  Lectures,  Department  of  Mathematics,  Univ.  of  Utah,  November  1986. 

5.  Awarded  Faculty  Fellowship,  1986-87,  University  of  Colorado. 

6.  Invited  Lecture,  Nonlinear  Workshop,  BYU,  Provo,  Utah,  March  1987. 

7.  Invited  Lecture,  Dept,  of  Mathematics,  BYU,  Provo,  Utah,  March  1987. 

8.  Invited  Lecture,  Dept,  of  Mathematics,  Utah  State,  March  1987. 

9.  Talk,  NATO  Advanced  Research  Workshop  on  Mathematical  Modeling  in  Combustion, 
Lyon,  France,  April  1987. 

10.  Invited  Lecturer,  Institut  Mathematique,  Universite  Catholique  de  Louvain,  Belgium, 
April-May,  1987. 

11.  Keynote  Lecturer,  Belgian  Mathematical  Society,  Hasselt,  Belgium,  May  1987. 

12.  Invited  Lecture,  Departmento  de  Matematicas,  Universidad  Complutense  de  Madrid, 
Spain,  May  1987. 


13.  Invited  Lecture,  Joint  USA/Japan  Workshop  on  Mathematical  Combustion,  Juneau, 
Alaska,  August  1987. 

14.  Invited  Talk,  Mathematics  Dept.,  Colorado  State  Univ.,  Fort  Collins,  November  1987. 

15.  Invited  Lecture,  Dept,  of  Mathematics,  Univ.  of  Pittsburgh,  Pittsburgh,  March  1987. 

16.  Invited  Lecturer,  Special  Year  on  Reaction-Diffusion,  Dept,  of  Mathematics,  Heriot-Watt 
Univ.,  Edinburgh,  Scotland,  May  1988. 

17.  Invited  Lecture,  Dept,  of  Phys.  Chemical,  Univ.  Leeds,  Leeds,  England,  May  1988. 

Outreach  activities  of  Prof.  D.  R.  Kassoy  include: 

1.  “Mathematical  Methods  for  Combustion  Modeling”:  A  series  of  12  lectures.  Department 
of  Aeronautical  Engineering,  University  of  Nagoya,  Nagoya,  Japan,  6/10-20/85. 

2.  “Unified  Theory  of  Low  and  High  Speed  Flames,”  University  of  Tokyo,  Dept,  of  Reaction 
Chemistry,  6/24/85. 

3.  Saitoma  University,  Mechanical  Engineering  Department,  Tokyo,  Japan,  6/25/85. 

4.  “Unified  Theory  of  Low  and  High  Speel  Flames,”  and  “Unified  Chain-Branching  &  Thermal 
Explosion  Theory,”  Tsukuba  University,  Engineering  Mechanics  Institute,  Japan,  6/26/85. 

5.  “Detonation  Wave  Initiation  by  Rapid  Energy  Deposition  at  a  Confining  Boundary,”  10th 
International  Colloquium  on  Dynamics  of  Explosions  and  Reactive  Systems,  S/P/85. 

6.  “The  Status  of  Thermal  Explosion  Singularity  Modeling,”  International  Workshop  on 
Mathematics  in  Combustion,  8/12/85. 

7.  “Compression  of  a  Confined  Inert  Gas  by  a  Piston,”  American  Physical  Society,  Fluid 
Dynamics  Division,  Tucson,  Arizona,  11/24/85. 

8.  “Compression  of  an  Inert  Gas  in  a  Cylinder  by  a  Piston,”  University  of  East  Anglia, 
Norwich,  England,  11/29/85;  and  Cranfield  Institute  of  Technology,  Cranfield,  England, 
12/4/85. 

9.  “Detonation  Wave  Initiation  by  Rapid  Energy  Deposition  at  a  Confining  Boundary,”  4th 
Army  Conference  on  Applied  Mathematics  and  Computing,  Ithaca,  New  York,  May  28-30, 
1986. 

10.  “Detonation  Wave  Initiation,”  International  Workshop  on  Mathematics  in  Combustion, 
Garmisch-Parten-Kirchcn,  Germany,  August  9-11,  1986. 

11.  “Subsonic  High  Speed  Reaction  Zone  Structure  for  Variable  Mole  Chemistry,”  Western 
States  Combustion  Institute  Meeting,  Tucson,  Arizona,  October  27-29,  1986. 

12.  “Compression  of  a  Confined  Inert  Gas  by  a  Piston,”  Fluid  Dynamics  Division  of  the  Amer¬ 
ican  Physical  Society,  November  23-25,  1986. 

13.  “Structure  of  a  Planar  High  Speed  Reaction  Zone,”  Fluid  Dynamics  Divison  of  the  Amer¬ 
ican  Physical  Society,  November  23-25,  1986. 

14.  “Initiation  of  Gasdynamics  by  Intense  Power  Deposition,”  Computational  Mathematics 
Group,  C.U.-Denvcr,  2/17/87. 

15.  “Detonation  Wave  Initiation  by  Energy  Deposition  at  a  Boundary,”  invited  presentation, 
Minisymposium  I,  SIAM  Conference  on  Numerical  Combusion,  San  Francisco,  3/10/87. 

16.  “Modeling  the  Early  Phases  of  Planar  Detonation  Wave  Initiation,”  2nd  ASME-JSME 
Thermal  Engineering  Joint  Conference,  Honolulu,  3/24/87. 

17.  “High  Speed  Reactioan  Zone  Structure  for  Variable  Mole  Chemistry,”  Mathematical  Sci¬ 
ences  Institute,  Cornell  University,  Ithaca,  5/4/87. 

18.  “High  Speed  Reaction  Zone  Structure  for  Variable  Mole  Chemistry,”  11th  International 
Colloquium  on  Dynamics  of  Explosions  and  Reactive  Systems,  Warsaw,  Poland,  8/5/87. 
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19.  “Detonation  Wave  Initiation  and  Related  Ignition  Phenomena,”  8th  International  Work¬ 
shop  on  Mathematics  in  Combustion,  Warsaw,  Poland,  8/4/87. 

20.  “Studies  Related  to  Planar  Detonation  Initiation  and  Straucture”  and  “Formulation  of 
Diffusive  and  Nondiffusive  Thermal  Explosion  Theory,”  JOint  ”USA/Japan  Seminar  on 
Mathematical  Modeling  in  Combustion  Science,  Alaska,  8/81/87  and  8/19/87. 

21.  “Diffusive  and  Nondiffusive  Thermal  Explosions,”  Mathematics  Department,  University  of 
Colorado,  9/3/87. 

22.  “Detonation  Wave  Initiation  by  Power  Deposition:  Numerics  and  Asymptotics,”  SIAM 
35th  Anniversary  Meeting,  10/13/87. 

23.  “A  Unified  Formulation  for  Diffusive  and  Nondiffusive  Thermal  Explosions,  University  of 
California  at  San  Diego,  4/20/88;  and  University  of  Southern  California,  4/22/88. 


3.  Publications 

Summaries  of  each  of  the  published  and  submitted  articles  follow: 

J.  Bebernes,  D.  Eberly,  and  W.  Fulks,  “Solution  profiles  for  some  simple  combustion  models,” 
Nonlinear  Analysis:  Theory,  Methods  and  Applications,  10  (1986),  217-228. 

Using  phase  space  techniques,  the  solution  shapes  for  the  Gelfand  problem  ~Au  =  Seu 
and  the  perturbed  Gelfand  problem  -Ati  =  £exp  6  >  0,  e  >  0  are  analyzed.  Both 

of  these  models  play  a  fundamental  role  in  the  mathematical  theory  of  thermal  explosions 
for  finite  rigid  and  gaseous  systems.  For  rigid  systems  the  physical  processes  are  determined 
by  a  pointwise  balance  between  chemical  heat  addition  and  heat  loss  by  conduction.  During 
the  inductive  period,  with  a  duration  measured  by  the  conduction  time  scale  of  the  bounding 
container,  the  heat  released  by  the  chemical  reaction  is  redistributed  by  thermal  conduction. 
As  the  temperature  of  the  container  increases,  the  reaction  rate  grows  dramatically.  Eventually 
the  characteristic  time  for  heat  release  becomes  significantly  smaller  than  the  conduction  time 
in  a  well-defined  hot  spot  embedded  in  the  system.  Then  the  heat  released  is  used  almost 
entirely  to  increase  the  hot  spot  temperature.  Both  the  models  studied  detect  this  hot  spot 
development  in  a  very  precise  manner.  For  example,  for  a  ball  J3j  C  JR3,  the  minimal  solution 
of  the  Gelfand  problem  is  bell-shaped,  that  is,  exhibits  the  hot  spot,  for  a  range  of  parameter 
values  6,  6  <  6  <  6*.  When  0  <  6  <  6,  the  physically  significant  minimal  solution  is  concave 
down  while  all  other  solutions  are  bell-shaped.  Here  6 *  is  the  Frank-Kamenetskii  critical  value 
and  6  is  a  uniquely  determined  parameter  value. 
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J.  Bebernes.  “Solid  Fuel  Combustion-Some  Mathematical  Problems”,  Rocky  Mountain  J. 
Math.,  16  (1986),  417-433. 

The  initiation  of  a  combustion  process  involves  a  myriad  of  complex  physical  phenom¬ 
ena  which  are  fascinating  to  observe  and  challenging  to  describe  in  quantitative  terms.  In 
general  one  is  concerned  with  the  time-history  of  a  spatially  varying  process  occurring  in  a  de¬ 
formable  material  in  which  there  is  a  strong  interaction  between  chemical  heat  release,  diffusive 
effects  associated  with  the  transport  properties,  bulk  material  motion  as  well  as  several  types 
of  propagating  wave  phenomena.  Mathematical  models  capable  of  describing  these  combustion 
systems  incorporate  not  only  familiar  reaction-diffusion  effects  associated  with  rigid  materials, 
but  those  arising  from  material  compressibility  as  well.  For  a  combustible  gas,  the  complete 
reactive  Navier-Stokes  equations  are  required  to  describe  the  phenomena  involved. 

In  this  paper,  we  shall  focus  on  the  initiation  and  evolution  of  thermal  explosion  processes 
in  rigid  materials.  In  this  situation  the  physical  processes  are  determined  by  a  pointwise  balance 
between  chemical  heat  addition  and  heat  loss  by  conduction. 

J.  Bebernes  and  D.  Kassoy,  “Solution  profiles  and  thermal  runaway,”  Lecture  in  Applied  Math¬ 
ematics,  24  (1986),  217-228. 

This  is  an  expository  description  of  how  thermal  runaway  occurs  for  the  induction  period 
model  of  a  high  activation  energy  thermal  explosion  in  a  bounded  container.  Assuming  a  slab 
geometry,  the  temperature  perturbation  0(x,t)  (solution  of  (*)  0t  =  A0  +  See)  blows  up  at  a 
single  point  x  =  0  as  t  — *  T  and  6(x,t )  — ►  0e(x)  for  0  <  x  <  1.  In  principle  0e(x)  is  found  from 
an  initial  value  numerical  solution  of  (*).  Assuming  0e(x)  as  a  known  quantity,  a  final  value 
theory  for  this  problem  is  described.  This  allows  one  to  predict  the  shape  of  the  temperature 
0(x,t)  in  a  neighborhood  of  the  singularity. 

J.  F.  Clarke,  D.  R.  Kassoy  and  N.  Riley,  “On  the  direct  initiation  of  a  plane  detonation  wave,” 
Proc.  Royal  Soc.  London  A408  (1986),  129-148. 

It  is  assumed  that  energy  is  transferred  at  a  rapid  rate  through  a  plane  wall  into  a  spatially 
uniform  and  initially  stagnant  combustible  gas  mixture.  This  action  generates  a  shock  wave,  just 
as  it  does  in  an  inert  mixture,  and  also  switches  on  a  significant  rate  of  chemical  reaction.  The 
Navier-Stokes  equations  for  plane  unsteady  flow  are  integrated  numerically  in  order  to  reveal 


the  subsequent  history  of  event.  Four  principal  time  domains  are  identified,  namely  ‘early’, 
‘transitional’,  ‘formation’,  and  ‘ZND/.  The  first  contains  a  conduction-dominated  explosion 
and  formation  of  a  shock  wave;  in  the  second  interval  the  shock  wave  is  reponsible  for  the 
acceleration  of  chemical  activity,  which  becomes  intense  during  the  ‘formation’  period.  Finally 
a  wave  whose  structure  is  in  essence  that  of  a  ZND  detonation  wave  emerges. 

J.  Debernes,  A.  Bressan,  and  D.  Eberly,  “An  asymptotic  description  of  blow-up,”  Indiana 
Univ.  J.  Math.  36  (1987),  295-305. 

The  purpose  of  this  paper  is  to  give  a  precise  description  of  the  asymptotic  behavior  of  a 
radially  symmetric  solution  u(x,t)  of 

(I)  ut  -  Au  =  eu 

in  a  neighborhood  of  the  blow-up  point  as  t  approaches  the  finite  blow-up  time  T  <  oo  provided 
x  €  Br  =  {x  :  |x|  <  J?}  C  Hln  and  n  =  1  or  2.  Giga  and  Kohn  recently  characterized  the 
asymptotic  behavior  of  the  solution  u(x,  t )  of 

{II)  ut  -  Ati  =  up ,  ficF 

near  a  blow-up  singularity  assuming  a  suitable  upper  bound  on  the  rate  of  blow-up,  provided 
n  =  l,2orn>3  and  p  <  .  For  fl  =  B\  C  IRn,  in  light  of  recent  a  priori  bounds  established 

by  Freedman  and  McLeod,  this  implies  that  the  solution  u{x,t )  of  (H)  with  suitable  initial  and 
boundary  conditions  satisfies 

{T  -t)pu{\x\,t)~+  P*3  as  t-+T~ 

provided  0  <  |z|  <  C(T  -  t)1/2  for  some  C  and  /?  =  ^y. 

For  (I)  we  prove  that  the  solution  u(x,t)  satisfies 

u(x,t)  -  In ■— - ♦  0 

v  ’  '  T-t 

uniformly  on  0  <  |x|  <  C(T  -  t )1/2,  C  >  0,  as  t  — *  T~ . 

Equation  (I)  is  the  ignition  period  model  for  the  thermal  explosion  of  a  nondeformable 
material  of  finite  extent  undergoing  a  single-step  exothermic  reaction.  This  model  neglects 
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the  consumption  of  fuel  and  describes  the  temperature  u  at  any  time  and  at  any  point  in  the 
bounded  container.  For  R  >  0  sufficiently  small,  the  temperature  u(x,<)  becomes  unbounded 
in  the  L 00 -sense  as  t  approaches  a  finite  blow-up  time  T  and  the  blow-up  occurs  at  a  single 
point,  the  center  of  the  ball  Br. 

Our  main  result  gives  a  precise  description  of  how  the  blow-up  asymptotically  behaves. 
The  proof  is  valid  only  in  dimensions  1  and  2.  The  real  difficulty  in  understanding  the  blow¬ 
up  for  (I)  lies  in  determining  the  nonincreasing  globally  Lipschitz  continuous  solutions  of  the 
associated  steady-state  equation 


on  [0, 00)  where  y'(0)  =  0,  y( 0)  =  a  >  0.  For  n  =  1  Bebernes  and  Troy  proved  that  the  only 
such  solution  is  y  =  0. 

Ivassoy-Poland  and  Kapila  in  earlier  papers  derived  (K)  from  an  asymptotic  final  time 
analysis  of  (I)  and  predicted  on  the  basis  of  numerical  calculations  the  existence  of  a  solution 
of  (Iv)  which  satisfies  the  asymptotic  condition  p(i?)  ~  -2lnr)  +  K0  as  7?  — ►  00.  This  implies  the 
existence  of  a  globally  Lipschitz  continuous  solution  of  (K).  Our  result  shows  that,  for  n  =  1  and 
2,  their  final  time  analysis  is  incomplete.  More  precisely,  we  remark  that:  i)  their  derivation 
of  the  asymptotic  condition  is  incorrect,  and  ii)  the  space  variable  in  the  inner  expansion 
rj  =  x/  \/T  —  t  is  stretched  too  much  so  that  no  information  concerning  the  spatial  behavior  of 
the  solution  is  retained  in  the  limit.  This  leaves  open  the  problem  of  finding  a  different  rescaling 
of  the  space  coordinates  for  which  in  the  limit  a  non-constant  behavior  of  solutions  can  still  be 
observed. 

J.  Bebernes  and  W.  Troy,  “Nonexistence  for  the  Kassoy  problem,”  SIAM  J.  of  Math.  Analysis , 
18  (1987),  1157-1162. 

When  blow-up  or  thermal  runaway  occurs  in  finite  time  for  the  thermal  explosion  of  a 
non  deformable  material,  the  developing  hot  spot  becomes  unbounded  at  a  single  point  of  the 
container  if  it  is  radially  symmetric.  These  supercritical  processes  which  are  characterized  by 
the  appearance  of  this  singularity  at  a  finite  time  T  have  been  considered  earlier  by  Kassoy- 
Poland  and  Kapila.  Using  computational  methods  for  symmetric  slab,  cylindrical,  and  spherical 
geometries,  they  predict  that  6(x,t )  becomes  unbounded  at  the  symmetry  point  x  —  0  at  time 


9 


T.  Elsewhere  6(x,t)  — *■  &e(x ),  x  €  fi,  x  ^  0,  as  t  — >  T.  By  applying  a  final- value  asymptotic 
analysis  at  T ,  they  predict  the  character  of  the  singularity  function  $i(x).  This  prediction 
is  based  on  the  existence  of  a  nontrivial  solution  of  y"  =  |y'  ey  -  1  =  0,  y'(0)  =  0,  and 
y(x)  ~  —2lnx  +  c  as  x  — *•  oc.  By  using  connectedness  and  shooting  technique  arguments,  we 
prove  that  no  such  solution  can  exist.  This  means  a  more  detailed  analysis  of  the  blow-up 
singularity  is  required  and  this  has  been  done  in  subsequent  papers. 

J.  Bebernes  and  W.  Troy,  “Existence  for  the  modified  Kassoy  problem,”  Proc.  Royal  Soc.  Ed¬ 
inburgh ,  106A  (1987),  131-136. 

In  this  paper,  the  differential  equation 

(1)  y"  ~  \y'  +  e*  -  1  =  0 

with 

(2)  y(0)  =  a,  y'(0)  =  f3 

is  considered.  Solutions  to  (l)-(2)  are  sought  which  have  the  asymptotic  property  y(x,  o,/3)  ~ 
—2lnx  +  Ka  as  x  -*  oc.  We  prove  that  there  exists  a  >  0  such  that  for  each  0  <  a  <  a,  there 
is  at  least  one  0(a)  <  0  such  that  (l)-(2)  has  a  solution  y(x,a,(3)  with  the  desired  asymptotic 
property. 

J.  Bebernes,  A.  Brcssan,  and  A.  Lacey,  “Total  blow-up  versus  single  point  blow-up,”  J.  Differ¬ 
ential  Equations,  73  (1988),  30-44. 

Traditional  thermal  explosion  theory  is  used  to  describe  reaction  initiation  in  condensed 
explosives  and  is  limited  formally  to  nondeformable  materials.  Kassoy  and  Poland  significantly 
extended  this  theory  to  develop  an  ignition  model  for  a  reactive  gas  in  a  bounded  container 
in  order  to  describe  the  induction  period.  During  this  induction  period  there  is  a  spatially 
homogeneous  pressure  rise  in  the  system  which  causes  a  compressive  heating  effect  in  the 
constant  volume  container.  Mathematically  this  compressibility  of  the  gas  is  expressed  by 
means  of  an  integral  term  in  the  induction  model  for  the  temperature  perturbation  0(x,t). 
This  model  is  given  by 

(D)  et-M  =  6e9  +  1—^  •  -Lr  /  9t(x,  t)dx 

7  voin  Jn 


Bebernes  and  Bressan  analyzed  earlier  ignition  model  (D)  for  a  compressible  gas  and 
proved:  For  any  6  >  0  and  any  7  >  1,  (D)  has  a  unique  classical  solution  6(x,t)  on  ft  X  [0,T) 
where  ft  C  Stn  is  a  bounded  container  and  T  <  +00  or  T  =  +  00.  In  the  latter  case,  6(x,t) 
blows  up  as  t  approaches  T.  If  6  >  i >fk >  the  Frank-Kamenetski  critical  value,  then  T  <  00  and 
blow-up  or  thermal  runaway  occurs  in  finite  time. 

In  this  paper,  a  description  of  where  blow-up  occurs  in  the  given  container  ft  is  given  for 
the  more  general  problem 

u(-Au  =  f(u)  +  g(t) 

( G )  u(  ,0)  =  4>(x),  x  €  ft 

u(x,i)  =  0,  x  £  3ft,  t  >  0. 

If  ft  is  a  bail  of  radius  R  in  lRn  and  blow-up  occurs  at  finite  time  T,  then: 

I)  If  /0T  g(T)dt  =  +00,  then  blow-up  occurs  everywhere. 

II)  If  /0r  g(t)dt  <  +00  and  f(u)  =  eu  or  /(u)  =  (w  +  A)p,  A  >  0,  p  >  1,  then  blow-up 
occurs  only  at  x  =  0. 

III)  If  </(0  =  ^  JQut(x,t)dx ,  1<  <  1,  and  /(u)  =  (u  +  A)p,  A  >  0,  1  <  p  <  1  +  2 /n, 
then  blow-up  occurs  everywhere. 

IV)  If  g(t)  =  /n  ut(x,t)dx,  K  <  1,  and  }{u)  =  eu,  then  blow-up  occurs  at  a  single 
point. 

V)  If  g(t)  =  fn  ut(x,t)dx  and  /(it)  =  (u  +  A)p,  A  >  0,  p  >  1  +  2/n,  then  blow-up 
occurs  only  at  1  =  0  provided  K  <  1  is  sufficiently  small. 

J.  Bebernes  and  D.  Ebcrly,  “A  description  of  self-similar  blow-up  for  dimensions  n  >  3,”  An- 
nalcs  de  I’Institut  Ilcnri  Poincare- Analyse  non  lineaire,  5  (1988),  1-21. 

The  purpose  of  this  paper  is  to  give  a  precise  description  of  the  asymptotic  behavior  for 
solutions  u(z,  t)  of 

(1)  ut  =  Au  +  f(u) 

which  blow-up  in  finite  positive  time  T.  We  assume  /(u)  =  up  (p  >  1)  or  /(u)  =  eu,  and 
z  C  Bn  =  {z  G  IRn  :  |^|  <  I?}  where  R  is  sufficiently  large  to  guarantee  blow-up. 
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Giga  and  Kohn  recently  characterized  the  asymptotic  behavior  of  solutions  u(z,t)  of  (1) 
with  f(u)  =  up  near  a  blow-up  singularity  assuming  a  suitable  upper  bound  on  the  rate  of 
blow-up  and  provided  n  =  1,2,  or  n  ^  3  and  p  5  ^±1.  For  c  jftn  usjng  recent  a  priori 
bounds  established  by  Friedman-McLeod,  this  implies  that  solutions  u(z,t)  of  (1)  with  suitable 
initial-boundary  conditions  satisfy 

(2)  (T-tfu(ztt)-+0fi  as  t-*T~ 

provided  \z\  <  C(T  —  tf)1/2  for  arbitrary  C  >  0  and  where  (3  = 

For  f(u )  =  du  and  n  =  1  or  2,  Bebernes,  Bressan,  and  Eberly  proved  that  solutions  u(z,t) 
of  (1)  satisfy 

(3)  u(z,  t)  4-  ln(T  —  t)  — ►  0  as  t  —*  T~ 

provided  \z\  <  C(T  —  t)1/2  for  arbitrary  C  ^  0. 

The  real  remaining  difficulty  in  understanding  how  the  single  point  blow-up  occurs  for  (1) 
rests  on  determining  the  nonincreasing  globally  Lipschitz  continuous  solutions  of  an  associated 
steady-state  equation 

(4)  y"  +  -  0  y'  +  F(y )  =  0,  0  <  x  <  oo 

where  F(y)  —  yp  —  (3y  or  ey  —  1  for  /(y)  =  yp  or  ev  respectively  and  where  y(0)  >  0  and 

y'(o)  =  o. 

For  F(y)  =  yp  -  (3y  in  the  cases  n  =  1,2,  or  n  ^  3  and  p  ^  ^7,  we  give  a  new  proof  of 
a  special  case  of  a  known  result  that  the  only  such  positive  solution  of  (4)  is  y(i)  =  f)0 .  For 
F(y)  =  ey  -  1  and  n  =  1,  Bebernes  and  Troy  proved  that  the  only  such  solution  is  y(z)  =  0. 
Eberly  gave  a  much  simpler  porof  showing  y(z)  =  0  is  the  only  solution  for  the  same  nonlinearity 
valid  for  n  =  1  and  2. 

For  3  <  n  5  9,  Troy  and  Eberly  proved  that  (4)  has  infinitely  many  nonincrcasing  globally 
Lipschitz  continuous  solutions  on  [0,oo)  for  F(y)  =  ev  -  1.  Troy  proved  a  similar  multiplicity 
result  for  (4)  with  F(y)  =  yp  -  fiy  for  3  <  n  5  9  and  p  >  s±l. 

This  multiple  existence  of  solutions  complicates  the  stability  analysis  required  to  precisely 
describe  the  evolution  of  the  time-dependent  solutions  u(z,  t)  of  (1)  near  the  blow-up  singularity. 


In  this  paper  we  extend  the  results  of  Giga-Kohn  and  Bebernes-Bressan-Eberly  to  the 
dimensions  n  t  3  by  proving  that,  in  spite  of  the  multiple  existence  of  solutions  of  (4),  the 
asymptotic  formulas  (2)  and  (3)  remain  the  same  as  in  dimensions  1  and  2.  The  key  to  un¬ 
raveling  these  problems  is  a  precise  understanding  of  the  behavior  of  the  nonconstant  solutions 
relative  to  a  singular  solution  of  (4)  given  by 

(5)  Se(i)  =  In 
for  f(u)  =  eu  and  n't  3,  and 

(6)  Sp(*)={-4/?[/?+±(2-n)]/:r2} 

for  /(u)  =  up  and  (3  -f  \{2  -  n)  <  0,  n  t  3.  This  will  be  accomplished  by  counting  how  many 
times  the  graphs  of  a  nonconstant  self-similar  solution  crosses  that  of  the  singular  solution. 

M.  A.  Birkan,  C.  K.  Law  and  D.  R.  Kassoy,  “Transient  decomposition-recombination  dynam¬ 
ics  of  dissociating  and  ioninzing  gases,”  Proc.  Royal  Soc.  London  A418  (1988),  331-343. 

A  generic  gas  AB,  initially  at  a  relatively  high  temperature,  undergoes  a  spatially  homoge¬ 
neous,  constant- volume  decomposition-recombination  reaction  represented  by  AB  +  M  A  + 
B  +  M.  The  complete  time-history  of  this  variable  temperature  reaction  is  calculated  by  em¬ 
ploying  high  activation  energy  asymptotic  analysis  developed  originally  for  thermal  explosion 
problems.  The  evolution  of  the  reaction  is  described  in  terms  of  an  initiation  period  with  small 
changes  in  the  mixture  temperature  and  composition,  a  longer  major  decomposition  period  dur¬ 
ing  which  most  of  the  conversion  of  AB  to  A  and  B  occurs,  and  an  extended  final  period  during 
which  recombination  becomes  important  as  the  system  relaxes  to  the  equilibrium  state.  Ex¬ 
plicit  timescales  are  derived  for  each  of  the  distinct  reaction  processes.  The  analytically  derived 
solution  agrees  quantitatively  with  the  numerical  solution  and  qualitatively  with  experimental 
results. 

D.  R.  Kassoy,  A.  K.  Kapila  and  D.  S.  Stewart,  “A  unified  formulation  for  diffusive  and  nondif- 
fusive  thermal  explosion  theory,”  Combustion  Science  and  Technology,  in  press  (1988). 

A  mathematical  model  is  developed  for  the  induction  period  of  a  parcel  of  slightly  warm 
reactive  gas  mixture  embedded  in  a  cooler  gaseous  environment.  High  activation  energy  asymp- 
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totic  methods  are  used  to  formulate  the  problem.  The  important  parameters  in  the  nondimen- 
sional  equations  are  ratios  of  characteristic  reaction,  acoustic  and  conduction  times  in  the 
thermally  distrubed  parcel  of  dimension  t1.  The  results  of  the  formulation  imply  that  a  spec¬ 
ified  physicochemical  system,  with  a  known  characteristic  chemical  time,  will  experience  a 
traditional  diffusively  dominated  thermal  explosion  if  the  dimension  V  is  sufficiently  small.  In 
larger  systems  transport  effects  are  negligible  and  the  induction  period  process  is  dominated 
by  reactive  gasdynamic  equations.  The  study  provides  a  unified  formulation  for  the  induction 
period  of  all  thermal  reaction  problems  and  specifies  the  relationship  between  classical  thermal 
explosion  theory  for  rigid  materials  and  compressible  gases  and  more  contemporary  efforts  to 
study  nondiffusive  thermal  reactions. 

D.  R.  Kassoy,  N.  Riley,  J.  Bebernes  and  A.  Bressan,  “The  confined  nondiffusive  thermal  explo¬ 
sion  with  spatially  homogeneous  pressure  variation,”  Combustion  Science  and  Technology ,  in 
press  (1988). 

A  solution  is  developed  for  a  nondiffusive  thermal  explosion  in  a  reactive  gas  confined  to 
a  bounded  container  fl  with  a  characteristic  length  l'.  The  process  evolves  with  a  spatially 
homogeneous  time-dependent  pressure  field  because  the  characteristic  reaction  time  t'R  is  large 
compared  to  the  acoustic  time  I'/Cq  where  Cq  is  the  initial  sound  speed.  Exact  solutions,  in 
terms  of  a  numerical  quadrature  are  obtained  for  the  induction  period  temperature,  density,  and 
pressure  perturbations  as  well  as  for  the  induced  velocity  field.  Traditional  single-point  ther¬ 
mal  runaway  singularities  are  found  for  temperature  and  density  when  the  initial  temperature 
disturbance  has  a  single  point  maximum.  In  contrast,  if  the  initial  maximum  is  spread  over  a 
finite  subdomain  of  17,  then  the  thermal  runaway  occurs  everywhere.  Asymptotic  expansions  of 
the  exact  solutions  are  used  to  provide  a  complete  understanding  of  the  singularities.  The  per¬ 
turbation  temperature  and  density  singularities  have  the  familiar  logarithmic  form  -  ln(t'  -  f') 
as  the  explosion  time  t'e  is  approached.  The  spatially  homogeneous  pressure  is  bounded  for 
single-point  explosions  but  is  logarithmically  singular  when  global  runaway  occurs.  Compres¬ 
sion  heating  associated  with  the  unbounded  perturbation  pressure  rise  is  the  physical  source  of 
the  global  thermal  runaway. 
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K.  Kirkkopru  and  D.  R.  Kassoy,  “Induction  zone  structure  of  a  high-speed  deflagration  with 
variable  mole  chemistry,”  The  Physics  of  Fluids,  in  press  (1988). 

The  induction  zone  characteristics  of  a  planar  subsonic  high-speed  reactive  flow  downstream 
of  a  specific  origin  is  investigated  theoretically  for  the  global  irreversible  reaction  F + Ox  -*  i>P. 
The  equation  of  state  for  the  reacting  gas  mixture  is  more  general  than  that  for  a  constant 
molecular  weight  gas.  Perturbation  methods  based  on  the  limit  of  high  activation  energy  are 
used  to  construct  the  general  parameter  dependent  analytical  solutions.  The  dependence  of  the 
ignition  delay  distance  on  the  kinetic,  stoichiometric  and  flow  parameters  is  discussed  in  detail. 
Significantly,  it  is  shown  that  the  maximum  ignition  delay  distance  when  the  chemical  heat 
addition  and  the  origin  values  of  parameters  are  fixed.  The  physics  and  length  scales  found 
from  the  perturbation  analysis  are  used  as  a  guide  in  generating  supporting  numerical  solutions. 

J.  Bebernes  and  D.  Eberly,  “Characterization  of  blow-up  for  a  semilinear  heat  equation  with  a 
convection  term,”  Quarterly  J.  Mech.  Math.,  accepted  (1988). 

Let  fi  C  JR"  be  a  bounded  domain  with  smooth  boundary  dtl.  Consider  the  initial¬ 
boundary  value  problem 

u<  =  Au  +  f(x,t,u),  (i,()£flx(0,T), 

u(x,  0)  =  «o(x)  >  0,  i6fl,  (1) 

u(x,t)  =  0,  (x,t)£dSl 

where  uo(Oil)  =  0.  Assuming  that  there  is  local  existence  of  a  classical  solution  of  (1),  a  recent 
active  area  of  interest  has  been  to  determine  and  characterize  how  solutions  of  (1)  cease  to 
exist.  For  a  large  class  of  nonlinearities,  solutions  cease  to  exist  by  becoming  unbounded  in 
some  norm  as  a  maximal  time  T  is  approached.  This  phenomenon  is  called  blow-up. 

In  particular,  the  following  problems  arise  naturally  when  considering  the  blow-up  of  so¬ 
lutions  and  have  been  considered  by  various  authors  for  special  nonlinearities  /. 

Problem  1.  Find  the  maximal  T  E  (0,oo]  for  which  a  classical  solution  to  (1)  exists  and 
determine  when 

a.  T  =  oo  and  u(x,t)  is  finite  onflx  [0,  oo],  or, 

b.  T  <  oo  and  ||  u  ||— »  oo  as  t  — *  T~  for  some  norm. 
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Problem  2.  Describe  the  set  of  blow-up  points  B  C  fi.  A  point  xq  G  fi  is  a  blow-up  point  for 
(1)  if  there  is  a  sequence  {(rm,im)}“=1  such  that  tm  — ►  T~,  xm  -» x0,  and  u(xm,tm)  — ►  oo  as 


Problem  3.  Describe  the  asymptotic  behavior  of  u(x,t)  as  t  — *  T~ .  In  particular,  at  a  blow-up 
point  xq,  describe  how  the  blow-up  singularity  evolves  with  the  ultimate  goal  being  to  describe 
the  space  profile  of  u(x,T )  in  the  neighborhood  of  an  isolated  blow-up  point. 

Problem  4.  For  t  >  T,  one  can  determine  what  happens  to  the  solution  u(x,t)  of  (1)  in  some 
“weak  solution”  sense. 

If  one  considers  the  more  general  semilinear  parabolic  problem 

ut  =  Au  +  F(x,t,u,  Vu),  (x,t)  6  D  X  (0,T) 


u(x,0)  =  u0(x)  >  0, 


u(x,  t)  =  o,  (x,  t)  e  on  x  (o,  r), 

then  Problems  1  through  4  remain  essentially  open.  For  the  special  case  where  F(x,  t,  u,  Vu)  = 
|u|p_1u  -  IVul’  with  p  >  1  and  q  >  1,  Chipot  and  Weissler  have  given  conditions  under  which 
solutions  blow  up  in  the  Loo-norm  in  finite  time  T.  In  this  case,  the  gradient  term  has  a  damping 
effect  working  against  blow-up,  so  these  results  become  more  delicate  to  establish.  A  critical 
splitting  of  q  given  by  gives  a  change  in  the  behavior  for  this  problem.  Chipot-YVeissler 
study  the  case  1  <  q  < 

In  this  note  we  consider  a  second  special  case  of  (2),  namely 

utAu  +  e“  -  |Vu|2,  (x,t)  G  D  X  (0,T), 

u(x,0)  =  u0(x)  >  0,  x  e  D,  (3) 

u(x,t)  =  0,  (x,t)  G  dfl  X  (0,T). 

which  can  be  viewed  as  the  limiting  case  of  the  critical  splitting  as  p  — ►  oo  in  the  problem 
considered  by  Chipot  and  Weissler.  Becaue  of  the  special  form  of  the  nonlinearity,  we  can  give 
very  precise  answers  to  the  stated  problems  1  through  4. 

The  keys  to  our  analysis  are  two  simple  observations.  First,  the  change  of  variable  v  =  1— e“ 
transforms  IBVP  (3)  into  the  equivalent  problem 

vt  =  Au  + 1,  (x,t)  e  n  x  (o,r), 

v(x,0)  =  u0(x),  x  G  D,  (4) 

v(x,  t)  =  0,  (x,t)  e3fix  (0,T), 


(4) 
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where  vo(du)  =  0  and  v0(i)  =  1  —  exp(tio(z))  <  1.  Second,  a  point  io  6  is  a  blow-up  point 
for  (3)  with  blow-up  time  T  if  and  only  if  v(x0,T)  =  1  at  some  first  time  T.  The  associated 
steady-state  problem  for  IBVP  (4)  is 


Ap  +  1  =  0,  x  €  SI, 
p(x)  =  0,  x  6  8SI. 


(5) 


It  is  well-known  that  IBVP  (4)  and  BVP  (5)  have  unique  solutions  t>(x,  t)  and  p(x ),  respectively. 
It  is  also  known  that  lim^oo  v(x,  t)  =  p(z)  uniformly  for  x  €  Cl. 

We  first  look  at  the  problem  of  when  blow-up  occurs.  As  we  have  observed,  blow-up  occurs 
at  xo  G  ft  at  time  T  for  IBVP  (3)  if  and  only  if  v(xo,T)  =  1  at  some  first  time  T.  Thus,  we 
need  only  determine  conditions  which  force  the  solution  t>(x,f)  of  (4)  to  attain  the  value  1. 
If  we  can  show  that  the  steady-state  solution  p(x)  of  BVP  (5)  attains  the  value  one  for  some 
x  €  ft,  then  blow-up  must  occur  for  IBVP  (3). 


K-  Kirkkopru  and  D.  R.  Kassoy,  “High-speel  reaction  zone  structure  for  variable  mole  chem¬ 
istry,”  SIAM  J.  Applied  Math.,  submitted  (1988). 

The  structure  of  a  high-speed  deflagration  downstream  of  a  specific  origin  is  investigated 
theoretically  for  the  global  irreversible  Arrhenius  type  reaction  F  +  Oz  — *■  vP.  Attention  is 
focused  on  the  effect  of  the  stoichiometric  coefficient  value  v  on  the  high-speed  reaction  zone 
characteristics.  Perturbation  and  numerical  methods  are  used  to  find  solutions.  In  the  rapid 
reaction  region  beyond  the  induction  zone,  there  is  a  strong  interaction  between  large  chemical 
heat  release  and  flow  compressibility.  When  the  chemical  reaction  causes  a  mole  reduction 
iy  <  2)  and  there  is  sufficiently  large  heat  release,  the  flow  velocity  reaches  a  maximum  and  then 
declines  while  the  temperature  increases  monotonically  throughout  the  process.  Significnatly, 
it  is  shown  that  the  flow  cannot  evolve  to  the  Chapman-Jouguet  (C.J.)  state  where  the  final 
local  Mach  number  is  unity  and  the  reactant  concentration  is  zero.  When  the  mole  number  of 
the  evolving  flow  is  constant  or  increasing  ( v  >  2),  the  flow  velocity  always  increases  and  a  final 
C.J.  state  can  develop  when  the  right  amount  of  chemical  heat  addition  is  available  for  a  given 
input  Mach  number  at  the  origin. 
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